CAS-KITPC/ITP-158 



Casimir Energy, Holographic Dark Energy and Electromagnetic 
Metamaterial Mimicking de Sitter Space 

Miao LiQ Rong-Xin MiaoQ and Yi Paiigp 

Kavli Institute for Theoretical Physics, 
Key Laboratory of Frontiers in Theoretical Physics, 
Institute of Theoretical Physics, Chinese Academy of Sciences, 
Beijing 100190, People's Republic of China, 
Interdisciplinary Center for Theoretical Study, 
University of Science and Technology of China, Hefei, Anhui 230026, China 



Abstract 

We compute the Casimir energy of the photon field in a static de Sitter space and find it to be 
proportional to the size of the horizon, the same form of the holographic dark energy. We suggest 
to make metamaterials to mimic static de Sitter space in laboratory and measure the predicted 
Casimir energy. 
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I. INTRODUCTION 



The dark energy problem is one of the central unsolved problems in fundamental physics 
and cosmology. Many phenomenological models have been proposed since the discovery 
of the accelerated expanding universe jlj, including the so-called holographic dark energy 
model Q. 

The value of the observed dark energy is comparable to the critical energy density which 
in the Planck unit is proportional to the Hubble constant squared. Motivated by this fact, 
the holographic dark energy model presumes that the energy density is inversely proportional 
to the square of the size of some cosmological horizon, and the choice of the horizon, for the 
model to work, is proposed in [2j] to be the event horizon. In a pure de Sitter space in which 
dark energy is a constant, one can choose a static coordinate system, thus the total energy 
of dark energy can be defined. The total energy in the Planck unit is proportional to the 
size of the horizon, unlike the usual result for the Casimir energy in the finite cavity in a flat 
spacetime, the latter is inversely proportional to the size of the cavity. Nevertheless, there 
has been no attempt to computing the Casimir energy in a static de Sitter space seriously, 
to the best of our knowledge. 

In this paper, we take the assumption that the UV divergent zero-point energy is somehow 
regularized to vanish due to some unknown mechanism, therefore the observed dark energy 
is solely due to the residual zero-point energy which is known to be the Casimir energy. 
However, we are faced with the challenge that why the Casimir energy is not inversely 
proportional to the size of the cavity. 

Motivated by the recent development in electromagnetic metamaterials, we will first try 
to convert the de Sitter space to a cavity of optical metamaterial, this is doable at least for 
the photon field. We shall find that the cavity is drastically different from a usual one in 
that the permittivity and permeability parameters both have divergent components tangent 
to the boundary, indeed these components are divergent near the boundary which is derived 
from the horizon of the de Sitter space. This fact encourages us to guess that the Casimir 
energy as a function of the size of the cavity behaves differently from that in a normal cavity. 
We carry out the calculation for the photon field and find that indeed after regularization 
there is still a divergent piece proportional to the size of the cavity, with the divergence 
cut-off by a UV scale. This UV scale in the metamaterials can be some intrinsic microscopic 
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scale, and becomes the Planck scale when transformed back to the de Sitter space. Thus, 
we find a piece of theoretical evidence for the holographic dark energy model. We propose 
to do some laboratory experiment to verify our result about the Casimir energy. This kind 
of experiment is interesting in two aspects. First, it will mimic cosmology, second, it will 
detect a unusual Casimir energy inversely proportional to some microscopic cut-off. 

There are two different designs of metamaterials, the first is based on the static de Sitter 
space with the usual coordinate system, the second is based on the static de Sitter where 
the radial coordinate is replaced by the proper distance. The Casimir energy in the second 
metamaterial assumes the precisely same form as in the de Sitter space, which is proportional 
to the size of the material thus is enhanced to be easily detectable. 

In the next section, we present the Maxwell equations in a curved spacetime. We discuss 
how to transform the de Sitter space into a cavity of metamaterial in sect. 3 and compute 
the Casimir energy in sect. 4, and discuss its form in a metamaterial filled cavity in sect. 5. 
We conclude in sect. 6. 

II. MAXWELL EQUATION IN STATIC CURVED SPACE TIME 

We start with Maxwell equations in a curved spacetime 

d^F vX] = 0, d^ v = 0, (1) 

where the first equation is the Bianchi identity satisfied by field strength F^ u , and H uv is 
defined by 

H** = y^F^ = y/^gg^F^. (2) 
The electric field Ei and magnetic field Hi are related to F^ u as 

E^F oi , H^-'-fH^, (3) 

In above equations, the Levi-Civita symbol e y is +1 for all even permutations of e 123 . It 
should be noted that Ei and Hi are spatial vectors. 

Especially, when the space time is static, we can express Eq. (TjQ) in a spatial covariant 
form with respect to optical metric 7 U defined by 

I 13 = -W j , (4) 
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where 7 U is the usually called optic metric. 

At this stage, Maxwell equations Eq.([T]) can be rewritten as 



ViE? = 0, ViH* = 0, (5) 
dtE* - ^LdjH k = 0, d t H i + ^d j E k = 0, (6) 



where V, denotes the covariant derivative with respect to 7 Jjf , and all the indices are raised 
and lowered by 7* J . 

Eq.([6]) can be rewritten in the following compact form 













: 









i% I _\=L\ __ I , (7) 

where the operator L is given by 

jkl 



If we define 

/ E 
H 

the operator L is Hermitian in terms of the following inner product 



(9) 



< 0i, 02 >= J ^ 3 xv^0|0 2 = J d 3 x(E u e ij E 2j + HufjFHv), (10) 

where 

e« = ^' = VtA 7 = det( 7fi ). (11) 

We will see that e and have a physical interpretation in terms of metamaterials in the next 
section. 

To make L Hermitian under the inner product (10) (which is the energy functional as it 
should be for 0i = 2 ), we need to specify boundary conditions if there is a boundary. We 
find the following boundary conditions 

y fy(E 1 xH 2 -H 1 xE 2 ) i m\ BM = (12) 

where rij is the normal vector of boundary and (E x H) % = -j=e l: * k EjH k . Then the hermiticity 
of L is guaranteed 

< 01, L0 2 > = < L0 X , 02 > • (13) 



By acting d t and L again on the two sides of Eq.(7) respectively, we obtain the equations 
of electric and magnetic field 

- d 2 E % = (V'V, - V 2 5j)E„ 

-d 2 H t = (V'V, - V 2 ^)F i5 (14) 
we read off the operator corresponding to L 2 from the above equations 

/ y^'Vi - V 2 5{ \ 

D = L = 1 . 15 

y V' V, - V 2 8l ) 

An eigenvalue of D is a nonnegative real number due to the Hermiticity of L. 
With Gauss law, we get two sets of equations taking the same form 

- d 2 V t = DjVj = - V 2 5l)V h W% = 0, (16) 

where Vi can be or H^. 

III. MIMICKING DE SITTER WITH METAMATERIALS 

Recently, the analogy between static curved space and inhomogeneous medium has been 
studied in [3|4j|. It is proved in [7j] that in geometric optic limit, light propagates along the 
same path in electromagnetic material and in static curved spacetime as long as the optic 
metric 7^ is related to permittivity and permeability by Eq. (1111) . 

One can see the analogy more clearly as follows. To solve Maxwell equation in inhomo- 
geneous medium conveniently, the often used quantities are electric displacement field D l 
and magnetic induction field B l defined as 

D i = e ij E 3 , B i = fi ij H 3 , (17) 

where £ u and /i u are permittivity and permeability specifying the electromagnetic properties 
of the medium. In terms of the newly defined quantities D l and H\ Maxwell equations in 
inhomogeneous medium and Cartesian coordinate system are 

dtD* = 0, 80 = 0, (18) 

dtD 1 — e ijk djH k = 0, d t B l + e ljk djE k = 0. (19) 



Since under spatial coordinate transformation e lJ and ^ are tensor density with weight -1, 
we can define two tensors using and ^ by 

£ ij -■■ /;'' 

y» = — , 7« = -Jl . (20) 

1 det(e*i)' ^ det(^) 1 ; 

For electromaganetic material with equal £ u and 7 U = 7 U . In this case, rewriting 
Eqs.( fT8l) and (fl9l) in terms of 7* J , ^ and ifj, we find that Maxwell equations in the metame- 
terial takes the same form as Eq.Q and Eq.Q if 7 U is identified with the optic metric 

rJ-3 _ 

By this analogy, one can design electromagnetic metamaterial corresponding a given met- 
ric. Along this line, theoretically, a new material called electromagnetic cloak was conceived 
3|-[6j. All parallel bundles of incident rays are bent around some region covered by electro- 
magnetic cloak and recombined in precisely the same direction as they entered the medium; 
Experimentally, a new class of electromagnetic materials are invented 8| which can be de- 
signed to have properties mimicking novel gravitational effects in Nature. Thus it is now 
conceivable that a material can be constructed whose permittivity and permeability values 
may be designed to vary arbitrarily. 

With the static de Sitter space 

ds 2 = -(1 - ^)d 2 t + (1 - r ^Y X d\ + r 2 d 2 n (21) 
The permittivity and permeability are 

sin 9 



e 



- = !T = r 2 sin 0, e ee = = ^ = ^ = - * . - , (22) 



where (r, 9, ip) stands for spherical coordinate. In terms of Cartesian coordinate, the per- 
mittivity and permeability are 

= i^TVZ^-V 5 - (23) 

The metamaterial can always be redesigned by making a coordinate mapping {3, 5|. By 
making the following mapping 

r — > r = L arcsin(r/L). (24) 

Then the permittivity and permeability becomes 

e" = ^ = L 2 ^P±\ sin 9, = ^ = -^- y = ^ = 1 . (25) 

cos(r/L) cos{r/L) cos{r/L) sin 9 



where (r, 6, ip) denotes the spherical coordinate. In terms of the Cartesian coordinate 



These two kinds of metamaterial mimicking the same de Sitter space, since they are 
related by a spatial coordinate transformation. However, the different permittivity and 
permeability will lead to distinct physical phenomenon in laboratory due to their different 
physical composition. 

The event horizon at r = L or f = ttL/2 now becomes the boundary of a cavity of 
metamaterial. As we shall see shortly, to metamaterial with e and \i as in Eq. (1231) . the leading 
term in the Casimir energy is a constant inversely proportional to certain microscopic length 
scale; while to metamaterial with e and fi as in Eq. (1261) . the leading term in the Casimir 
energy in this cavity is proportional to L, we encourage experimenters to design such cavity 
to measure the Casimir energy, the result will have important implication for dark energy. 

IV. CASIMIR ENERGY OF ELECTROMAGNETIC FIELD IN STATIC DE SIT- 
TER SPACE 

In Maxwell theory, the electric sector and magnetic sector contribute the same amount 
to the Casimir energy, as the equations governing them are identical. In the following, we 
will focus on the electric sector. 

Consider the solution of electric field taking the form 



where the absolute value of u should be interpretated as the energy with respect to time 
t. Upon quantization, the eigenvalue u contributes to the zero-point energy. The Casimir 
energy of electric field is the infinite sum 



Substitutting Eq. (f2T|) into Eq.( TT6|) we obtain the eigen equation of operator D with 
eigenvalue u 2 




(26) 




(27) 




(28) 



D\Ej = u 2 Ei. 



(29) 
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Then in the framework of zeta function regularization, the problem of computing the Casimir 
energy is converted into the problem of computingthe heat kernel of D, the Casimir Energy 
can be extracted from the heat kernel as follows [9J: 

-^Casimir = lim ^{£ reg (+e) + £ reg (-e)}, (30) 

e— »0 Z 



where 



^eg(6) = f C(-~ + e )> ( 31 ) 



with the zeta function associated with D defined by 

,2s roo 

/ dt t s -Hi\e- tD ). (32) 
Jo 



C(s) 



In above expressions, a scale "//" is introduced to keep the zeta function dimensionless; while 
the prime in tr' indicates that zero eigenvalues are not included. In coordinate representa- 
tion, (x\e~ tD \x) is called the heat kernel denoted by K(t,x,x) satisfying the heat equation 
related to D as 

d t K(t,x,x) + DK{t,x,x) = 0. (33) 

To solve this equation, we need to specify appropriate boundary conditions. 

At this moment, with above preparation, we can calculate the Casimir energy of electro- 
magnetic field in static de Sitter space time. This space time is described by the following 
metric 

2 2 

ds 2 = -(1 - r —)dH + (1 - ^r'dV + r 2 d 2 n (34) 

where L is the de Sitter radius. However, the effective metric relevant to our calculation is 
the optic metric 7 U appearing in the Maxwell equation with the following form 

2 2 

ds\ = lt] dx l dx^ = (1 - yy\(l - JjT 1 ^ + r2rf2fi )> (35) 

where the lower index "3" is used to emphasize that this metric is three dimensional. An 
interesting observation is that 7^ actually describes an anti de Sitter space. One can see 
this clearly by performing the following transformation on r coordinates: 

r = . r (36) 
^/l + r 2 /L 2 

In terms of f, 

~2 

ds\ = (1 + ^)~ l d 2 r + f 2 d 2 Q. (37) 



In the case of static de Sitter space time, the boundary condition Eq. ffl2|) for retaining the 
Hermiticity of operator L in Eq.f U5p becomes 



(E w H 2( t > — H w E 2( j ) )\ r= L-8 = 0, (38) 

where the boundary is chosen at r = L — 5 and 5 is a cut off to avoid the divergence of energy 
and its meaning will become clear later (see later for an explicit calculation of energy). Thus 
to keep Hermiticity, we choose 

E e \ r=L _ s = E lf \ r=L _ s = (39) 

as the boundary condition. It is similar to the case of electromagnetic field in a spherical 
conductor shell. Therefore in the following, we will concentrate on the problem of finding 
the heat kernel corresponding to D defined with respect to 7^ in Eq.( )35|) and restricted to 
the space spanned by its eigenfunctions satisfying the Gauss law and boundary conditions. 

V*Ei = 0, E e \r=L-6 = E v \r= L -s = 0. (40) 

Usually, in a curved manifold the heat equation Eq. (1331) is hard to solve completely. However, 
the trace of the heat kernel has the following asymptotic expansion for small t, especially in 
three dimensions as 

tr'(e- tD ) = (-i-p{ f a/ + f b n t n ) + o(t N ) }, (41) 

where a n are bulk terms independent of boundary condition, composed by polynomial of 
Riemann tensor, Ricci tensor, Ricci scalar, and their covariant derivatives; b n are boundary 
terms comprised of the curvature, second fundamental form and their derivatives on the 
boundary. The sum runs over half-integers, but a n vanishes for half-odd-integers. Above 
asymptotic expansion implies that ((s) has a pole structure 9|: 

where f(s) is an entire analytic function of s, but in general, one has little information about 



it {9]. In the previous version of this work, we calculate the Casimir energy with above 
formulas. It is found that the leading term comes from a surface integral on the boundary 
horizon. Based on this expression, one may think that the Casimir energy is localized near 



the boundary horizon, but this may be a misunderstanding. For instance, in the frame work 
of heat kernel expansion, the Casimir energy between two infinite conductor plane in flat 
space will turn up as a boundary term depending on the permittivity and permeability of the 



conductor appearing to be localized on the boundary However, L. Ford's result [10J shows 
that in this case the Casimir energy distributes homogeneously between the two conductor 
planes. 

To avoid this misunderstanding, We will seak for a compact expression for the heat kernel. 
It is recalled that the optical metric describes an Euclidean AdS3 space and the heat kernel 
defined with respect to it can be solved by utilizing Harmonic analysis 11J. To vector field 
the corresponding heat kernel satisfying the transverse condition and vanishing at the AdS3 
boundary was given by [l2] (In our convention, the AdS3 boundary corresponds to r = L 
or f = oo). From Eq. (HOI) , it is obvious that the result of [3] amounts to the 5 — > 



12| as a good 



limit of the heat kernel in our case. In other words, we can use the result of 
approximation to the heat kernel sought by us up to the order S/L. Based on that result 
we write down the contracted coincident limit of the heat kernel as 



7 ii ^ i (x,x,t) 



t 



(47Tt)l L 2 ' 

With above result, tr'(e _tD ) can be evaluate by 



o(S/L). 



(43) 



tr'(e" tD ) = j y/yflKyfaxrfcPx 

= -L\- s ~ ln(^)) x ^ (1 + 2^)e-^ 2 + o(l), (44) 

(47Tt)2 Is 

where o(l) denotes terms <C k. Substituting this result into Eqs. fl32|) and fl30|) . we obtain 

q i t t or 

£ca Simir = Y^O^ 2 " 7 - ^)){- 5 - 1 Mj)) + o{l/L), (45) 

where 7 is Euler constant and ip(— |) = r'(— |)/r(— |). In above result, 5 is a coordinate 
cut off so the physical meaning is not clear. To express it in terms of the physical cut off, 
we recall that in static de Sitter space the physical distance between r = L — 5 and r = L 
is given by 

A/ phy = L aicsm{r/L)\ r r = L L _ s . (46) 
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By requiring A/ P h y to be the Planck length l p , we derive 



5 = L(l - cos(^)) 



L 7/ 2L 2 ' 

In terms of the physical cut off, the Casimir energy of electromagnetic field is 



(47) 



Eca 



in 



;(ln^ 2 - 7 - *(r\))i^ - i M^)) + o(l/L), 



(48) 



where the dominant term proportional to L/ti takes the same form as in the holographic 



dark energy model 



2]. Compared with previous result given in 
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14j, it is found that 



the leading term in the vacuum energy density of static de Sitter space is proportional to 



,2\-2 



which can also generate a leading term proportional to L/fi after integrated over 



the bulk of static de Sitter space. Finally we note that the running renormalization In /i 2 
can be absorbed into ti or be determined experimentally 13]. 



V. CASIMIR ENERGY OF ELECTROMAGNETIC FIELD IN ELECTROMAG- 
NETIC METAMATERIAL MIMICKING DE SITTER 

We have designed two kinds of metamaterial mimicking de Sitter in sect. 3. For an elec- 
tromagnetic metamaterial with e and fi taking the value of Eq. (l23j) in laboratory, coordinate 
r has been the physical distance between origin and r, since metric measuring the physical 
distance is flat. In this case, the physical IR cut off is chosen to be L — d, and the Casimir 
energy is given by 

3 1 1 1 2T 

£cas imir = —(\n^ 2 - 7 - ^— 2 ))(- d - 1 H-j)) + o(l/L), (49) 

where d is the counterpart of the Planck length in metamaterial, its value depends on the 
detailed chemical components and structure of the material which should be determined by 
experiment. We note that when L ^> d, 

3 11 

^Casimir- — (hi /i 2 -7 -^(-^))^- (50) 
l07T 2 d 

L is usuaully much larger than d, the leading term is much larger than the usual term 
proportional to 1/L (for example, if d is nanometer, and L is 1cm, then the leading term is 
about 10 6 times of the second term, the usual term which is almost undetectable). 
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Alternatively, to metamaterial with e and \x taking the form as in Eq. (l26|) . the Casimir 
energy is 

3 /, 2 ,/ 1 w r-i / sin(f/L) 1 .1 + sin(f/L),,. 
So— = -(lV-7-*(- 5 ))L '(J^lWi^lllw^V^W 

(51) 

where f stands for physical distance between origin and f in this coordinate system and the 
size of the cavity is tcL/2. So distinguished from the first case, ttL/2 — d is the physical IR 
cut off. d has the same meaning as in previous case. Usually L^> d 

^Casimir « 7^0^ ~ 7 ~ V>(-^))4- (52) 

87r 2 or 

It is remarkable that in this metamaterial, the Casimir energy has the same form as in 

its gravity analogy, except that some microscopic scale d takes the place of Planck scale. 

VI. CONCLUSION 

It is a surprising result that the Casimir energy in a de Sitter space is not inversely 
proportional to the size of the horizon as for a usual finite cavity, rather it is proportional to 
the size of the horizon. However, if one is to expect that the Casimir energy is dark energy 
or at least a part of dark energy, this is a desired result, the normal answer in a cavity is 
just too small to be relevant in cosmology. 

Although we only computed the Casimir energy due to the photon field, it is expected 
that other fields will have the same form of Casimir energy. Since we suggest to construct 
a cavity of metamaterial to do experiment in laboratory, it is enough to know the form of 
the photon field Casimir energy, since this is what to be measured in laboratory. However, 
because the corresponding Casimir force is combined with the elastic force, it seems difficult 
to measure the predicted Casimir energy immediately. 

Aside from the issue of dark energy, to use metamaterials to mimic cosmology by itself 
is an exciting future direction, we look forward to developments in the future. 
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